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ABSTRACT 


The manner in which the steady flow divides at an arterial 
junction was experimentally established in terms of the significant 
dimensionless parameters: the entrance Reynolds number, Re, the ratio, 
8B, of theveree branch diameter to the main branch diameter and the 
non-dimensional angle of branching, a. 

The type of bifurcation examined was the "side branch", 
where a single branch leaves the parent trunk at an angle dependent 
on the anatomical position. 

The ratio, y, of the flow rate in the side branch to the 
flow rate in the main branch was found to increase with a decrease 
of the Reynolds number and the angle of branching and with an increase 
of the ratio of diameters, 8. 

A visualization technique proved the existence of two inter- 
dependent separation regions, one in each branch, responsible for the 
variation of the mass flow ratio as function of Re,8, 8. The forma- 
tion, growth and shedding of vortices in the separated region of the 
main branch and the double-helicoidal flow in the side branch have 


been observed. 
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6 angle of branching 
r = U/u', ratio of maximum to mean longitudinal velocity at 


the bifurcation entrance 


v fluid kinematic viscosity 

fe) fluid density 

oO =4Le/Re, dimensionless entrance length 

Tq thickness of the separated region in the side branch 


(see Figure 4.22) 
thickness of the separated region in the main branch 


(see Figure 4.22) 
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CHAPTER I 
INTRODUCTION 


1.1 Introductory Remarks 

As the human knowledge increases, one tries to apply it to 
a better understanding of oneself. This is particularly true in 
medical and paramedical research for which curing diseases and saving 
lives are the two main goals. 

Since diseases and defects of the human cardiovascular 
system remain an important cause of death, researchers are expending 
considerable energy to understand this complex system. The flow of 
blood through the vessels is one of the aspects of this system well 
suited for the engineer to investigate, due to the fact that it is 
essentially a problem in fluid mechanics, although a rather complicated 
one. 

Because of the high degree of complexity of the circulatory 
system, a systematic application of the physical laws and equations 
valid in viscous flow would be a considerable task, indeed. However, 
the use of experimental data as background and justification of theories 
based on some simplifications of the physiological model, can lead to 
important results and will contribute to the extension of human know- 


ledge. 
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2 “Literature Survey | 

One way of studying a highly complicated problem is to 
simplify it as far as possible, to study the new problem and then add 
more and more complexity until the study of the real problem can be 
done. For that reason, the blood flow study is very much indebted 
to the French physician Poiseuille who, for his blood flow study, 
turned to the study of fully-developed steady flows of Newtonian 
fluids in straight rigid tubes. 

The fully-developed flows appear when the effects of ie 
ends of the tube are negligible. H.L. Langhaar! and H.B. Atabek® 
considered these end effects in developing an inlet length theory 
for steady and periodic flows in straight rigid tubes, respectively. 
ett. Gerrard? experimentally investigated pulsating turbulent flow in 
a tube. 

N.R. Kuchar and S. Ostrach* studied entrance effects associ- 
ated with laminar flows of viscous fluids in circular elastic tubes. 
Later, in 1966, B.H. Ruterbories and S. Ostrach? measured velocity 
profiles and pressure distributions of pulsating flow in a flexible 
tube. 

The problem becomes even more complex in the case of branched 
vessels and, consequently, most of the work in this area is done with 
the help of experimental models. 

The study of steady flows through branches is of interest 
not only to the blood flow research but also in air-conditioning or 


air flows through the lungs. This depends on the type of bifurcation 
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and the kind of results required. For example, A. Vazsony1® gave 
formulas permitting calculation of losses in duct branches, S.F. 
Gilman’ studied the pressure losses of dividing flow fittings and 
T.J. Pedley, R.C. Schroter and M.F. Sudlow® studied the flow and 

the pressure drop in systems of repeatedly branching tubes. They 
used Y-type bifurcations as experimental representation of the lungs 
airways, obtained velocity profiles and energy dissipations through 
these bifurcations. 

When one wants to study pulsatile flow through a bifurca- 
tion, a new problem appears due to pulse wave reflections. One of 
the earlier attempts to analyse this phenomenon was done by R.L. Evans 
et Ae who compared reflections in a mode] with those in an animal. 
In 1966 J.D. Martin and M.E. clark!® experimentally and theoretically 
analysed wave reflections in branched flexible conduits. R.H. Cox!! 
and G.S. Malindzak!@ made more sophisticated theoretical and in-vivo 
analyses of the wave reflections, respectively. 

The influence of the branching on the critical Reynolds 
number has been determined by Lid: Krovetz!?, who found a lower critical 
Reynolds number in branched tubes than in straight tubes. The hydro- 
dynamic stability appeared independent of the varying angles of branching. 

The occurrence and the severity of the separated Flow, as 
a function of the flow rate, have been studied by D.J. Schneck and 
W.H. Gutstein!*, They noted a direct correlation between the occurrence 
of artherosclerotic plaques and the angle of bending and branching. 


The correlation between turbulent flow, thrombosis and arterial lesions 
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has been also pointed out by J.R. Mitchell and C.J. Schwartz !>. 


A more complete investigation of pulsatile flow in a side- 
16 
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branch bifurcation has been done by C.M. Rodkiewicz and D.H. Howell 
who systematically studied the manner in which the flow divides at an 
arterial bifurcation as a function of the significant dimensionless 
parameters: Reynolds number, unsteadiness parameter and velocity 
fluctuation parameter. They found that for certain values of these 
parameters, more flow goes into the side branch. This leads to the 
important conclusion, that the heart, activated by an impulse from the 
brain, can supply more blood to the preferential area of tne human 


body, by changing the values of these three parameters. 


1.3 Objectives of the Thesis 

In order to have a better understanding of the flow through 
a side-branch bifurcation and find why and when the mass flow ratio 
y is greater than one, it has been decided, from the results obtained 


by C.M. Rodkiewicz and D.H. Howel1!® 


, to study the steady-flow charac- 
teristics in an arterial junction. The parameters for this study will 
be: 

- the entrance Reynolds atta 

- the geometry of the bifurcation, 

- the smoothness of the inside edges of the bifurcation. 

Qualitative as well as quantitative results will be produced 

giving the variation of y, the occurrence and the severity of the 


separated flow and the location of the separation point, as functions 


of the total flow rate. 
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CHAPTER II 
DERIVATION OF THE PARAMETERS FROM A 
MATHEMATICAL APPROACH 


2.1 Introductory Remarks 


As previously explained, this thesis is mainly concerned 
with the in-vitro study of the blood flow characteristics in arterial 
junctions. 

The present study is limited to the steady case of a homo- 
geneous incompressible Newtonian fluid of constant viscosity flowing 


through a side-branch bifurcation. 


2.2 Governing Equations 


2.2.1 The General Continuity Equation is: 


in which op is the fluid density and q the fluid velocity at a point. 
Applying the conditions of steadiness (9/dt = 0) and incompressibility 


(9 = constant), the equation reduces to: 


V-q=0 (2.2) 
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2.2.2 The General Navier-Stokes Equations for a fluid with constant 


viscosity can be expressed as follows: 
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D L “ 
Dt i NDE ai vW(Veq)tvV 


i 
my 


q (2.3) 


where F and p represent the body force and the static pressure at the 
point under consideration and vy is the fluid kinematic viscosity. 
Writing these equations for the steady flow of an incompressible 
fluid, we finally obtain, using equation (2.2) and neglecting the 
body force (F = 0): 


vp t+ v v°q. (2.4) 


2.2.3 The Energy Equation can be written, for the sake of complete- 


ness of this chapter, in its general form as follows: 


V(kVT) + vo (2.5) 


in which e is the specific internal energy, q the heat generated per 
unit mass, k the fluid thermal conductivity, @ the dissipation function 
per unit mass and T the local temperature. 

In the case of the steady flow of an incompressible fluid 
with constant properties and without heat generation, equation (2.5) 


becomes: 
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2.3 Derivation of the Parameters 
2.3.1 Introductory Remarks 
The purpose of this section is to derive the parameters 
describing the flow entering the bifurcation; indeed, the flow through 
the bifurcation depends on these parameters. This necessitates further 
simplifications of the physiological model: 
1. Upstream of the bifurcation, the vessel is considered to be 
Straight and of constant circular cross-section and of rigid walls. 
2. The system is assumed to be dynamically and geometrically 
axisymmetric. 
2.3.2 Reduction of the Governing Equations 
The cylindrical polar coordinate system, as indicated in 
Figure 2.1, is the best one to be used in the case of a flow through a 
straight pipe and can aiso be used to describe the flow through the 
bifurcation. 
Written in this coordinate system, the governing equations 
(2.2), (2.4) and (2.6) become, taking into consideration the axial 
symmetry upstream of the bifurcation: 


1. Continuity equation 
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2. Navier-Stokes equations 
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3. Energy Equation 
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2.3.3 Deduction of the Parameters by Non-Dimensionalization Technique 


The coordinates, the velocity components, the local temper- 
ature and pressure are non-dimensionalized by use of the following 


transformations: 


Ne agp Decal 
r = Rr ve VV Ket) 
See p = AP p 


where the barred quantities are non-dimensional; U is the maximum 
longitudinal velocity component, V is the reference velocity in the 
radial direction, AP and AT are the maximum pressure drop and temperature 


difference in the entrance region respectively, R is the tube radius 
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and L the approximate length of the entrance region. The introduction 


of these quantities into the equations (2.7) through (2.10) allows the 


following set to be written: 
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In order that both terms of the non-dimensionalized con- 


tinuity equation (2.12) are of the same order of magnitude, set the 
factor [VL/UR] equal to one; that is: 


pak 
a: (2.16) 
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Substituting this relationship into the previous equations yields the 


following: 


he ea (2.17) 
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From equations (2.18), (2.19) and (2.20) the following set 


of non-dimensional parameters can be deduced: 
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2. Reynolds Number: Re = S 
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3. Euler Number: Eu = ai 
oU 
(220) 
4. Peclet Number: Pe = cek 
ye 
5. Eckert Number: E = CAT 


For the purpose of this thesis, which is not concerned with 
the pressure or temperature distribution throughout the bifurcation, 
the flow at the bifurcation inlet is fully described by only two of 


these dimensionless parameters, that is: 


le the Lengtm Ratio; Le =F 


2. the Reynolds Number: Re - 8 
The Reynolds number characterizing the flow through a bifurcation is 
calculated at the entrance to the bifurcation. If the entry length 
upstream of the bifurcation is sufficient, the velocity profile at 
the bifurcation inlet is parabolic. This allows us to express the 
Reynolds number in terms of the average longitudinal velocity com- 


ponent u'. The average velocity is equal to half of the maximum 


velocity: 
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and the Reynolds number becomes 


where D is the tube diameter. 


(2623) 
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CHAPTER III 
APPARATUS AND INSTRUMENTATION 


3.1 Requirements 


The purpose of the experimental model is to enable the 

study of the flow characteristics in arterial junctions. As explained 
in Chapter I, the present thesis confines itself to the steady case of 
blood flowing through a side-branch bifurcation. In order to deter- 
mine the range of study, it is necessary to summarize the blood 
characteristics related to the present work. The following data has 
been collected by C.M. Rodkiewicz and D.H. Howe11!®; 

1. Blood density varies from 1.048 to 1.066 g/cm. 

2. Blood relative viscosity with respect to distilled water 
varies from 3.5 to 5.4. 

3. Blood average velocity through a 2.1] cm diameter aorta is 
scacn/S. 

4. Human aorta diameter varies from ].5 to 2.1] cm. 
Substitution of this data into equation (2.23) yields the following 


entrance Reynolds number range: 
1300 < Re < 2040 (3.1) 


The bifurcation geometry, as illustrated in Figure 3.1, 
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Figure 3.1 


(a) with sharp inside edaes 


(b) with rounded inside edges 


Experimental Bifurcations 
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leads to the introduction of two further dimensionless geometrical 
parameters: oa related to the angle of branching, 6, and 8 as the ratio 
of the side-branch internal diameter, d, to the mainline-branch internal 


diameter, D; that is: 


a = 999 (322) 
g= 4 | (3.3) 
The range of study for o and 8 are: 
Fsa <] (3.4) 
0.40<8<1 (3.5) 


The main objectives of the experiments are: 

1. To study the dependence of y, ratio of the side-branch to 
the mainline-branch flow rates, on the three dimensionless numbers Re, 
a and B and to establish when and why y is greater than one. 

2. To establish, for each value of the couple (a,8), the Reynolds 
number at which the boundary layer starts to separate from the wall and 
then to observe the variation in position of the separation point 
when the Reynolds number varies. 


3. To determine the influence of the smoothness of the bifurca- 


tion edges on the flow. 
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4. To draw some conclusions regarding future research on the 
blood flow. 
Finally, in order to have a better understanding of the 
flow evolutions, the experimental range for the Reynolds number has 


been extended to: 


500 < Re < 5000 (3.6) 


3.2 Description of the Apparatus 
gezat ‘Choice of. the Experimental Fluid 


Taking into consideration the simplifications of the physio- 
logical model, already used as basis of the mathematical approach, 
the flow of blood through a given bifurcation, for a given Reynolds 
number at its entrance, should be similar to the flow of any incom- 
pressible Newtonian fluid through a similar bifurcation, for the same 
entrance Reynolds number. For this reason, as well as for ease of 
utilization, water has been chosen as the experimental fluid. 

In fact, in order to decrease the surface tension of the 
water with respect to the acrylic building material of the bifurcations, 
a wetting agent* and an anti foam** have been added to the water. 

Such a solution has a viscosity slightly different than water. The 
*Photo-Flo 200 Kodak, Edmonton Photo Supply, Edmonton, Alberta. 


**Anti-Foam Kodak, Edmonton Photo Supply, Edmonton, Alberta. 
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variation of its kinematic viscosity with the temperature has been 
measured with a Cannon-Fenske viscometer and reported in Appendix A. 
3.2.2 Overall Description of the Apparatus 

Figures 3.2 and 3.3 are the layout and a photograph of the 
apparatus, respectively. The apparatus can be divided into three 
different sections: 

1. The flow generation unit 

2. The flow measurement unit 

S.umine- test section 
These three sections are supported by a dexion frame and connected to- 
gether by copper pipes or flexible tubings of 1.25 inch internal dia- 
meter (this diameter was chosen for ease of observation). 
3.2.3 Flow Generation Unit (Figure 3.4) 

The required steady flow is obtained by using a centrifugal 
pump* in series with a damping tank and a Filter**. The damping tank 
is made of stainless steel and has an internal diameter of 10 inches 
and a height of 30 inches. Five sixths of this tank is filled with 
water, so that the remaining air cushion absorbs the possible fluctua- 
tions of the pump output. The outlet of the tank is located 20 inches 
above the bottom and is connected to the filter. 

The total flow rate and consequently the entrance Reynolds 


ee a 


*Motor-Mount Centrifugal Pump Deming, Unit No. 137, Electrical In- 
dustries Ltd. , Edmonton, Alberta. 


**Honeycomb Filter Tube No. EI5RIOA, B. Guthrie Engineering Company 
Ltd. , Edmonton, Alberta. 
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number is controlled by a 1-1/2" gate valve (A) at the outlet of the 
pump and by a 3/4" seat valve (B) located on the by-pass (valves A 
and B are shown on Figure 3.2). This flow generation unit is capable 
of delivering a maximum of 40 U.S. gallons per minute. 

3.2.4 Flow Rate Measurement Unit (Figure 3.5) 

As shown on the apparatus layout (Figure 3.2) and on the 
photographs (Figure 3.5), this unit has been divided into two dif- 
ferent sections: 

1. The slow flow measurement tanks. 

2. the flowrator meters.* 
Low flow rates are determined, in each branch, by simultaneously 
closing the drains of the slow flow measurement tanks and measuring 
the elapsed time to collect a known quantity of water (varying from 
1 to 5 liters). These tanks, made of clear plastic, were calibrated 
with extreme care. 

When the flow is fast enough, the branch flow rates are 
measured by the flowrator meters, which have been specially cali- 
brated to have an accuracy of one percent of the maximum flow. Each 
flow meter can measure flow rates up to 1.52 U.S. gallons per minute. 
3.2.5 Test Section (Figure 3.6 

In order to enable observation and photography, the bifurca- 


tions, shown in Figure 3.8, were made of acrylic. Two kinds of bi- 


*Flowrator Meter Fisher and Porter, Model 10A3565S Calibrated, Fisher 
& Porter (Canada) Ltd., Edmonton, Alberta. 
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Figure 3.5 Flow Rate Measurement Unit 


(b) Flowrator Meters 
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Figure 3.7 Generation Unit for Hydrogen Bubbles 
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Figure 3.8 Two Typical Bifurcations 


Figure 3.9 Disassembled Bifurca 


aS 

furcations have been built: 

1. For the first kind, with sharp inside edges, the two branches 
were drilled in a 1-1/2" thick block of transparent plastic. 

2. For the second kind, with rounded inside edges, two separate 
half-bifurcations were milled in 3/4" thick acrylic sheets (see 
Figure 3.9). The inside edges were then rounded by hand with curva- 
ture radius of 3/16" to 3/8" as illustrated in Figure 3.1. The two 
half bifurcations, positioned by two pins, were finally bolted 
together and sealed with chloroform. 

The use of plastic blocks and sheets gives a flat upper 
outside surface to the bifurcations. This offsets serious optical 
distortions which would have been created by a curved air-plastic 
interface (see Figure 3.10). Table 3.1 shows the geometries and the 
corresponding dimensionless geometrical parameters of the tested 
bifurcations. 

The inlet to the test section is a 1-1/4" internal diameter 
plastic pipe. Each bifurcation branch discharges into a constant 
head tank (see Figures 3.2 and 3.6). The two tanks have been designed 
to approximate the discharge of water from the bifurcation into two 
infinitely large reservoirs. They are identical, maintain the same 
water level at the two outlets of the bifurcation - 9" of water above 
the axis of each branch - and contain approximately 19 liters of water. 
The water overflows, through a 1 inch wide opening, from each reservoir 


to the corresponding measuring tank. 
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$72 <6) yEntrance Length 
The upstream branch is 8'7" long, which gives a length ratio 


equal to: 
Le = 164.8 (3/78) 


The coefficient o = 4 Le/Re, as defined by H.L. Langhaar| corres - 
ponding to the critical Reynolds number for a straight pipe (see 


H. Schlichting!/), is: 
o = 0.143 (3.8) 


According to reference 1 for this value of o, the ratio of the maxi- 
mum velocity to the mean velocity, at the bifurcation entrance, is 


given by: 
ay eee 


This means that the velocity profile, at the bifurcation inlet, approxi- 
mates the parabolic profile within 3.52. 
One can notice that 
1. This approximation shrinks down to 0.50% when the Reynolds 
number is 1650. 
2. In order to approximate the parabolic profile within 0.50% 


for a Reynolds number of 2300, the entry length would have to be 16'4", 
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which was not practical because of space limitations and would be of 


doubtful benefit. 


3.3 Flow Visualization 
One of the best qualitative as well as quantitative methods 
of visualizing water flow is the "Hydrogen Bubble Technique". 


3.3.1 Introductory Remarks on Electrolysis 


A typical electrolysis circuit is shown in Figure 3.11. 


Electrolytic Solution 


Figure 3.1] Electrolysis 
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When the battery or the DC power supply is operating, the electrons, 
drained away from the anode, are crowded onto the cathode. The 
electric circuit can only be closed if, simultaneously, a reduction 
and an oxidation occur respectively at the cathode and the anode. 
Using water as the electrolytic conductor, when a metal 
with a positive oxidation potential is chosen as a cathode material, 


the only possible reduction at this electrode will be: 
2H bubte os H,(gas) (3.10) 


The formation of oxygen gas from water by the equation (3.11) has an 
oxidation potential equal to -1.23V. If the anode is made of a metal 
with an oxidation potential greater than -1.23V, no oxygen can appear 


at the anode, since the oxidation of the anode is more likely to occur. 


2H,0 > 0,(gas) + 4H’ + de (3711) 


Consequently the use of a tungsten cathode and a copper anode wil] 
yield to the following reactions, since the oxidation potentials of 
tungsten and copper are + 0.12V and -0.34V, respectively: 


Anode: Cu+ 20H” + Cu(OH), + 2e (3.12) 


Cathode: PH canes Ho (gas) (3413) 
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Water Equilibrium: H,0 2 H’ + OH” (3.14) 


3.3.2 Hydrogen Bubble Technique 


The hydrogen bubble technique utilizes the electrolysis of 
water to introduce hydrogen into the flow. These bubbles, formed at 
the fine-wire cathode (0.002" diameter), are swept off by the flow. 
This 0.002" diameter wire produces 0.001" diameter bubbles !8, which | 
have a maximum vertical velocity, due to buoyancy effect (which in- 
creases with the size of the bubbles) of 0.001] fps. and reach a 
velocity of 98 percent of the fluid velocity within 1 ms, as explained 
by W. Davis and R.W. Fox !9, 

The wire has, due to its size, an insignificant influence 
on the flow itself. This constitutes one of the best advantages of 
this method, which also provides a great flexibility to yield qualita- 
tive and quantitative conclusions. Qualitative results can be ob- 
tained by continuously producing bubbles. This visualizes the pathlines, 
which are also streamlines and streaklines for a steady flow. If the 
bubbles are generated by electrical pulses, with a known frequency, 
one can simultaneously measure the velocities over a whole plane at a 
given instant. 

The success of this visualization method is highly dependent 
on the correct illumination. A sharply collimated sheet of powerful 
light is necessary as well as dark backgrounds and suppressicn of other 


light sources. 
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3.3.3 Experimental Equipment 

A special electrodes housing (Figure 3.12) has been built to 
hold both cathode (a 0.002" diameter tungsten wire) and anode (a copper 
ring). This housing has been designed so that: 

1. It perfectly matches the internal diameters of both up- 
stream pipe and bifurcation inlet. 

2. It can rotate around its axis in order to generate sheets 
of bubbles in any of the diametrical planes at the inlet to the bi- 
furcation. 

3. It enables the change of wire. 

The electrolysis DC current is brought to these electrodes 
from the generation unit (see Figures 3.7 and 3.13), which allows the 
operator to have complete and precise control of the frequency and the 
duration of the pulses producing the bubbles. The main part of this 
generation unit is the constant DC potential power supply*. Its 360V 
are applied between the two electrodes each time as long as the 
electronic switch (see diagram on Figure 3.14) is turned on by the 
pulse generator** driven by an oscillator’. The oscillator regulates 
the frequency of the pulses, while the pulse unit controls their dura- 
tion. Finally a reversing switch has been installed so that the wire 

*Lambda Regulated Power Supply, Model C 882M, Lambda Electronics Corp. 
**Unit Pulse Generator, Type No. 121/7C, General Radio Company. 


“Oscillator 2C-2MC, Type No. 1310A, General Radio Company. 
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36 
can be cleaned by reversing the electrodes polarity. 
Photographs of the streamlines are taken with a reflex 
camera* loaded with a high speed film** while the proper illumination 
is obtained by use of two slide projectors (total power 450W) giving 


a collimated sheet of light in the same plane as the bubbles sheet. 


‘Miranda Camera with Auto-Miranda Lens 1:1.9, f = 50 mm and Extension 
~Jube. 


**Kodak Tri-X Pan, ASA 400. 
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CHAPTER IV 
EXPERIMENTS, RESULTS AND DISCUSSION 


4.1 Experimental Procedure 
4.1.1 Preliminary Procedure 


This preliminary procedure, done after exchange of the filter- 
ing element is: 

1. If not already in place, the appropriate bifurcation is in- 
stalled. 

2. The pump is started, while valves A and C are opened and 
valve B is closed (see Figure 3.2). The damping tank is then under 
atmospheric pressure. 

3. When the damping tank is 5/6 full of water, the valve C is 
closed. 

4.1.2 Actual Experimental Procedure 

The two first steps of the following procedure are to be dis- 
regarded, when using the preliminary procedure. This actual experi- 
mental procedure is: 

1. The appropriate bifurcation is installed. 

2. The pump is started, while valvesB and C are closed and valve 
A is opened. 


3. When the water almost overflows the constant head tanks valve 


A is closed and the pump is stopped. 
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4. After the system reaches equilibrium, the levels of the two 
constant head tanks are adjusted to be exactly the same. This con- 
stitutes the most important step, as far as coherence in the results 
is concerned. 

5. The pump is then started and the flow is set as slow as 
possible by the use of the flow control valves A and B (see Figure 3.2). 

6. The electronics for hydrogen bubble generation are turned on. 

7. The flow rates in the two bifurcation branches are measured 
using the slow flow measurement tanks. 

8. The temperature of the solution is measured in one of the 
two constant head tanks. 

9. Observation of the flow through the bifurcation is made and 
photographs are eventually taken. 

10. The total flow rate is increased. 

11. Operations 7 through 10 are repeated until the flow is fast 
enough to be measured by the flow rator meters with sufficient precision. 

12. Operations 8 through 11 are repeated until the flow rates are 
no longer within the range of study. 

13. In order to drain the test section, valves A and B are closed 
and opened, respectively, while the pump is stopped. 

The entrance Reynolds number is calculated from the values of 

the total flow rate and of the temperature of the solution by the use 


of either coefficient a or coefficient b, given in Appendix A. 
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4.2 Experimental Results 

The experimental data are given in Appendix B. 
4.2.1 Mass Flow Ratio, y, as Function of Reynolds Number, Re 

Figures 4.1 through 4.4 and 4.5 through 4.8 give the variation 
of the mass flow ratio, y, with respect to the Reynolds number, Re, as 
a function of the diameter ratio, 8, for the sharp and round edged bi- 
furcations, respectively. 

Figures 4.9 through 4.12 and 4.13 through 4.16 give the vari- 
ation of y with respect to Re as a function of a for the sharp and 
round edged bifurcations, respectively. 

| Figure 4.17 gives a comparison between the two kinds of bi- 

furcations for a = 5/9. 
4.2.2 Mass Flow Ratio, y, as Function of the Bifurcation Geometry 

Figures 4.18 and 4.19 give the variation of y with respect to 
8B as a function of Re for the sharp and round edged bifurcations with 
B = 5/9 , respectively. 

Figures 4.20 and 4.21 give the variation of y with respect 
to a as a function of Re for the sharp and round edged bifurcations with 


B = 0.60, respectively. 
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Figure 4.1 Sharp Edged Bifurcation 


Mass Flow Ratio y versus Reynolds Number Re 


aS. aa Function, of: Beiormo-= 1/3 (6 <3 30°) 
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Figure 4.2 Sharp Edged Bifurcation 


Mass Flow Ratio y versus Reynolds Number Re 
as a Function of 8 for a = 5/9 (@ = 50°) 
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Figure 4.3 Sharp Edged Bifurcation 


Mass Flow Ratio y versus Reynolds Number Re 


as a Function of 8 for a = 7/9 (@ = 70°) 
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Figure 4.4 Sharp Edged Bifurcation 


Mass Flow Ratio y versus Reynolds Number Re 
as a Function of 8 for a = 1 (6 = 90°) 
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Figure 4.5 Round Edged Bifurcation 
Mass Flow Ratio y versus Reynolds Number Re 
as a Function of 8 fora= 1/3 (@ = 30°) 
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“Figure 4.6 Round Edged Bifurcation 


Mass Flow Ratio y versus Reynolds Number Re 


as a Function of 8 for a = 5/9 (6 = 50°) 
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Figure 4.7 Round Edged Bifurcation 


Mass Flow Ratio y versus Reynolds Number Re 


as a Function of 8 for a = 7/9 (6 = 70°) 
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Figure 4.8 Round Edged Bifurcation 
Mass Flow Ratio y versus Reynolds Number Re 
as a Function of 8 for a= 1 (6 = 90°) 


| ae 
re dhias 


y Ratio 


3000 4000 9000 


1000 2000 


Reynolds Number 


Figure 4.9 Sharp Edged Bifurcation 


Mass Flow Ratio y versus Reynolds Number Re 
as a Function of a for 6 = 1.00 
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Figure 4.10 Sharp Edged Bifurcation 
Mass Flow Ratio y versus Reynolds Number Re 
as a.Function of a. for 6 = 0:80 
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Figure 4.11 Sharp Edged Bifurcation 


Mass Flow Ratio y versus Reynolds Number Re 
asa FUNCEION 07-0) Ton tBe=20260 
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Figure 4.12 Sharp Edged Bifurcation 
Mass Flow Ratio y versus Reynolds Number Re 
as a Function of a for B = 0.40 
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Figure 4.13 Round Eaged Bifurcation 
Mass Flow Ratio y versus Reynolds Number Re 
as a Function of a for 8 = 1.00 
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Figure 4.14 Round Edged Bifurcation 
Mass Flow Ratio y versus Reynolds Number Re 
as a Function of a for 8 = 0.80 
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Figure 4.15 Round Edged Bifurcation 
Mass Flow Ratio y versus Reynolds Number Re 
as a Function of o for 6 =10700 
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Figure 4.16 Round Edged Bifurcation 
Mass Flow Ratio y versus Reynolds Number Re 
as a Function of a for g = 0.40 
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Figure 4.17 Comparison Between the 
Two Kinds of Bifurcations 


Mass Flow Ratio y versus Reynolds Number Re 
as a Function of 8 for a = 5/9 (@ = 50°) 
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Figure 4.18 Sharp Edged Bifurcation 
Mass Flow Ratio y versus Diameter Ratio B as a 
Function of Reynolds Number Re for 
a = 5/9 (6 = 50°) 
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Figure 4.19 Round Edged Bifurcation 
Mass Flow Ratio y versus Diameter Ratio 8 as a 
Function of Reynolds Number Re for 
oa = 5/9 (@ = 50°) 
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Figure 4.20 Sharp Edged Bifurcation 
Mass Flow Ratio y versus Non-Dimensional Angle oa 
as a Function of Reynolds Number for 8 = 0.60 
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Figure 4.2] Round Edged Bifurcation 


Mass Flow Ratio y versus Non-Dimensional Angle a 
as a Function of Reynolds Number for g = 0.60 
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4.2.3 Description of the Flow 
| Figure 4.22 represents the flow in a bifurcation. This figure 

was sketched from photographs and observations. 

Figures 4.23 and 4.24 show vortex formations in the separated 
region of the main branch and helicoidal flows in the side branch of a 
bifurcation, respectively. 

The severity of the separated regions, for a given bifurcation, 
as a function of the Reynolds number, is given in Figures 4.25 and 4.26. 

The dependence of the flow characteristics on the diameter 
ratio is shown in Figure 4.27, whereas their dependence on the branching- 


angle is studied in Figures 4.28 and 4.29. 


ane fF paar ine Pisa elt Wy + wort ont 


cehi@t tains ont hn 
6 to dooerd eke any eye WaT 1 sate ort fiat bn ors 
. 


Jnehiao all Pdopiaety. 98 we lt nota tease aint %6 ai 
yh Uk GS aarp cal _ Be od reo stem oi 


iva 4 (othe sae 8o Joga aL) ° f fe  2hAwiny m tem: oa) at 


Top view \e~ 


Section A-A 


Separated | 


separated region 


region 


Figure 4.22 Physical Description of the Flow 
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4.3 Observations and Discussions 
heel Mass Flow Ratio, y 

The study of the curves given in Figures 4.1] through 4.17 
leads to the following conclusions: 

1. The mass flow ratio, y, for a given bifurcation, is a de- 
creasing function of the Reynolds number. Since an increase of the 
Reynolds number increases the fluid momentum upstream of the bifurca- 
tion, less fluid will flow through the side branch. 

2. For a Reynolds number between 2000 and 3000, the value of 
the mass flow ratio, y, for a given bifurcation, depends on the flow 
history. Figure 4.10 shows that y is higher when the flow is increased 
from laminar to turbulent than when it is decreased from turbulent to 
laminar. This constitutes the transition region of the bifurcation. 

3. A given increase of the diameter ratio, 8, with everything 
else unchanged, will result in an increase of the mass flow ratio, y. 
This increase is a decreasing function of 8, as also shown in Figures. 
4.18 and 4.19. 

4. A given decrease of the angle of branching, 8, with every- 
thing else unchanged, will result in an increase of the mass flow ratio, 
y. This increase is a decreasing function of ala = @/90°) as also 
shown in Figure 4.20 and 4.21. 

5. The previous conclusions are valid whether the bifurcation 
has sharp or rounded inside edges. The general behaviour of the flow 
is then the same for the two kinds of bifurcations. However, for a 


Reynolds number less than 2000, the mass flow ratio, y, everything else 
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being unchanged, is slightly higher for a round edged than for a sharp 
edged bifurcation (see Figure 4.17). This, consequently, applies to 
the physiological range of Reynolds number (see paragraph 3.1) which 


is: 
1300 < Re < 2040 (4.1) 


This difference is probably due to the fact that the point of separation 
in the side branch is pushed downstream into the side branch for a 
round edged bifurcation with respect to the corresponding sharp edged 
bifurcation. 
4.3.2 Separation Regions 
Photographs and observations clearly prove the existence of 
two interdependent separation regions: one in the main branch and one 
in the side branch. If we defined their Heres by Tp and Tg: 
respectively, according to Figure 4.22, Figures 4.25 through 4.29 show 
that, whenever Tp increases, Tq decreases and vice versa. A growth of 
the separated region in the main branch will then restrict the flow 
in the main branch and will, consequently, increase y. 
The thickness, Tp» of the separated region in the main branch 
increases when: 
1. the Reynolds number decreases (Figures 4.25 and 4.26) 
9. the diameter ratio, 8, increases (Figure 4.27) 


3. the angle of branching decreases (Figures 4.28 and 4.29). 
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This is then the physical explanation of the Previous conclusions on 
the mass flow ratio, y. 

Inside the separated region in the main branch, back flows 
generate vortices, as shown in Figures 4.22, 4.23 and 4.24 (e and f). 
These vortices are not limited to the diametrical plane of the bifur- 
cation. In fact, they have been observed to present a banana-shape, 
which enables particles to be transferred from the main branch into 
the side branch, as sketched in Figure 4.22. 

Inside the side branch, the flow becomes double-helicoidal 
(see Figures 4.22 and 4.24). This kind of flow is typical of flows 
in curved pipes, as explained by H. Schlichting. |/ 
4.3.3 Limitations of the Apparatus 

The apparatus limits the accuracy of the results to approxi- 
mately + 2%. These results have been found to be reproducible within 
+ 2% for Reynolds numbers greater than 1000 (see Figure 4.5). This re- 
producibility is directly proportional to the care taken to obtain 
precisely the same level for the two constant head tanks, since the 
flow through the bifurcation is highly dependent on the resistance in 
each branch. 

For Reynolds numbers less than 1000 the reproducibility be- 
comes questionable, due to surface tension effects in the constant head 


tanks. Also, below that Reynolds number (Re = 1000), the visualization 


technique becomes useless, due to the growing importance of the buoyancy 


effect on the hydrogen bubbles. 
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Finally, from the observations and the photographs, no definitive 
conclusion can be drawn as far as the jocation of the separation point 


in the main branch is concerned. 
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CHAPTER Y 
CONCLUSTONS 


5.1 Summary 


In conclusion , for the physiological range of Reynolds number, 
the experimental findings can be summarized as follows: 

1. The mass flow ratio, y, is a decreasing function of the 
entrance Reynolds number. 

2. The mass flow ratio, y, iS an increasing function of the 
diameter ratio, 8. 

3. The mass flow ratio, y, is a decreasing function of the 
angle of branching, 06. 

4. The mass flow ratio, y, is slightly higher for a round edged 
bifurcation than for the corresponding sharp edged bifurcation, at 
the same entrance Reynolds number. 

5. There exist two interdependent separation regions, one in 
each branch. 

6. An increase of the mass flow ratio, y, is due to both an 
increase of the thickness, Tp» of the separated region in the main 
branch and a decrease of the thickness, Tq of the separated region 
in the side branch. 

7. Generation, growth and shedding of banana shaped vortices 


have been observed in the separated region of the main branch. In the 
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15 
Side branch, double-helicoidal flow has been observed. 
8. The flow through the bifurcation is highly dependent on the 


resistance downstream of the bifurcation. 


5.2 Suggestions for Further Work 


As previously explained, the apparatus has been found in- 
adequate to give valuable quantitative as well as qualitative informa- 
tion for slow flow (Re < 1000). Reducing the width of the Openings 
at the top of the two constant head tanks will eliminate the surface 
tension problems which appear at slow flow, and using the visualization 


20 for small fluid velocities will 


technique recommended by D.J. Baker 
be an excellent alternative for the hydrogen bubble technique. 

This visualization technique uses a pH indicator. If, on 
reverse position, (see Figure 3.13), the d.c. voltage applied between 
the two electrodes is kept sufficiently low, no hydrogen bubbles can 
be produced at the copper ring. Moreover, if the aqueous solution con- 
tains 0.01% of thymolsulphonephtalein, and is titrated to be on the acid 
(or yellow) side of the end-point (pH = 8.0), the production of 
hydroxy] jons (OH ) near the wire will change the colour from yellow to 
blue. Blue rows will then be generated in place of the previous rows 
of bubbles. Since the thymol blue molecule remains an ion in solution, 
no buoyancy effects are present. This method requires an illumination 


by a sodium arc lamp and cannot be used for flows faster than 2"/sec, 


since the coloured fluid is swept away from the wire more rapidly than 


it can be formed in visible amounts. 
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At very slow flow, levelling the two constant head tanks be- 
comes of primary importance due to the principle of connected vessels. 
Thus, before undergoing any slow flow tests, one has to find a means 
of being certain that the openings of the two constant head tanks are 
exactly at the same level. 

Not only the study of very slow flow should lead to interesting 
results but also the study of very fast flow should indicate whether y 
goes to zero when the Reynolds number increases, or not. The present 
apparatus can also be modified for this type of study by enlarging the 
Openings at the top of the two constant head tanks and using larger 
flowrator meters. 

For the range of Reynolds number (1000 < Re < 5000) for which 
the actual apparatus gives valuable results, four kinds of studies could 
be undertaken. 

1. The testing of other types of bifurcation such as Y-type or 
curved bifurcations. | 

2. Obtaining the velocity profiles throughout the bifurcation by 
measuring the distances separating two rows of bubbles directly on the 
photographs. This necessitates, in order to avoid any error due to 
optical distortions, to superimpose two photographs: one of the actual 
experiment and the other of a scale in the plane of the wire inside the 


bifurcation filled with water. 


3. The influence of the resistance downstream of the bifurcation 
on the flow through the bifurcation. 
4. The study of the frequency and pattern of vortex shedding, as 


a function of the Reynolds number and the bifurcation geometry. 
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APPENDIX A 
VISCOSITY OF THE EXPERIMENTAL FLUID 


The experimental fluid is water to which a wetting agent 
and an antifoam have been added. This solution has a viscosity slightly 
different than water. 

The kinematic viscosities of the solution and of the tap- 
water of the University of Alberta have been measured with a Cannon- 
Fenske viscometer.* Table A.1] and Figure A.1 give the variations of 
these viscosities with the temperature. 

It appears that the addition of the wetting agent and the 
antifoam affects the viscosity of the water by, at the most, 1.30%. 

When the temperature of the solution varies, its viscosity 
varies, so does the Reynolds number for a given flow rate. This 


Reynolds number can be expressed in terms of the total flow rate 


(Q) + Q)- Since 


OR eth ul Tr 7. CAS) 


the Reynolds number becomes 


*Cannon-Fenske-Ostwald Calibrated YViscometer No. 50 J657. 
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One can write this as 

Re = a (Q, + Q)) (A.3) 
when the flow rates are measured in liters per minute, or as 

Re = b (Q, + Q,) (A.4) 
when the measurements are made in U.S. gallons per minute. 


Table A.2 and Figures A.2 and A.3 give the variations of 


a and b with the temperature of the solution. 
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Table A.1 
Kinematic Viscosity (centistokes) as 


a Function of Temperature 


Temperature Tap-Water Solution 

2c oF 

20 68 1.015 1.002 

40 104 0.669 0.662 

60 140 0.487 0.482 

80 176 0.380 0.376 
Table A.2 


Variation of a and b with the Temperature 


eshte, 22 eee ee ee 


Temperature a b 

°C . (liters /mn)7| (U.S. gpm) 7! 
20 68 668 2536 

40 104 1012 3838 

60 140 1388 5267 

80 176 1781 B735 
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Figure A.1 Kinematic Viscosity of the 


Solution versus Temperature 
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Figure A.2 Value of a versus Temperature 
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Table B.1] 
Dependence of y on Re for 


a Sharp Edged Bifurcation with 


8 = 1.00 anda = 1/3 (6 = 30°) 
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Table B.2 
Dependence of y on Re for 


a Sharp Edged Bifurcation with 


B=1anda= 5/9 (@ = 50°) 
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Table B.3 
Dependence of y on Re for 
a Sharp Edged Bifurcation with 
8 = 1.00 anda= 7/9 (6 = 70°) 
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Table B.4 
Dependence of y on Re for 
a Sharp Edged Bifurcation with 
8B = 1 anda = 1 (@ = 90°) 
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Table B.5 
Dependence of y on Re for 
a Sharp Edged Bifurcation with 
8 = 0.80 anda=1/3 (6 = 30°) 
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Table B.6 
Dependence of y on Re for 
a Sharp Edged Bifurcation with 


B = 0.80 and a = 5/9 (8 = 50°) 
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Table B.7 


Dependence of y on Re for 


a Sharp Edged Bifurcation with 


8 = 0.80 anda = 7/9 (6 = 70°) 
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Table B.8 
Dependence of y on Re for 
a Sharp Edged Bifurcation with 
B= 0.80 anda=1 (6 = 90°) 


Los 
ie 
ik 
ibs 
Ie 
ih 
is 
ie 
Is 
i 
OF 
0. 
Ue 
OF 
UE 
0. 
OF 
0. 


Oe 
OF 
OF 
a 
0. 
03 
Q. 
0. 
Or 
0. 
0. 
We 
0. 
0. 
ae 
0. 
OR 


94 


t 
tl 
a 
_— 
) oF 
I 
i 
» A 
v A et ‘y 
' ihe 
Hi ft | Re: a% " wey 
Se. bisa " 
ih uf Pi 4 ‘ivi 
7 i » hie 7 “peor | 
i a : 


poe) ae 
ot ae ae 


Table B.9 
Dependence of y on Re for 
a Sharp Edged Bifurcation with 
p= 0500 anda = 1/3 (0 = 30°) 
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Table B.10 


Dependence of y on Re for 


a Sharp Edged Bifurcation with 


B = 0.60 and a = 5/9 (6 = 50°) 


2560 
2/750 
2930 
3130 
3320 
3600 
3780 
4000 
4100 
4190 
4300 
4490 
4740 
4950 


0.870 
0.87] 
0.878 
.856 
. 840 
2705) 


0 

0 

) 

0. 
OF 
Ue 
OF 
OF 
0. 
0. 


=) 


96 


a rennet 


“oo *: iy yt vi 


i 
i 7 | ie] y An ine A 
a F , : Ce ool ao * rd 
MN 4 Ar e4 a4 Te 
iy ra ; { i. aS 7 bi 


eet a y om ot 


= a ee = 
= Ss > ee wes 


P e ‘€ x 
= a % 


Table B.11 
Dependence of y on Re for 
a Sharp Edged Bifurcation with 
pa= 0 607andve = 7/916 =) /0°) 


2520 0. 


2750 
3000 
3200 
3450 
3700 
3840 
4000 
4200 
4320 
4380 
4500 
4600 
4800 
5000 


0 

0 

0 

0. 
OE 
0. 
UF 
Qi: 
0. 
OF 
0. 


845 


. 848 
.845 
Gc 
817 
cae 
./80 


97 


a \ 
pan 


Jail 


i i a 


Table B.12 
Dependence of y on Re for 
a Sharp Edged Bifurcation with 
B = 0.60 anda = 1 tole 90°) 
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Table B.13 
Dependence of y on Re for 
a Sharp Edged Bifurcation with 
8 = 0.40 and a = 1/3 (6 = 30°) 
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Table B.14 
Dependence of y on Re for 
a Sharp Edged Bifurcation with 
B = 0.40 and a = 5/9 (0 = 50°) 
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Table B.15 
Dependence of y on Re for 
a Sharp Edged Bifurcation with 
8 = 0.40 anda= 7/9 (6 = 70°) 
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Table B.16 
Dependence of y on Re for 
a Sharp Edged Bifurcation with 
B = 0.40 anda=1 #£2(6 = 90°) 
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Table B.17 
Dependence of y on Re for 
a Round Edged Bifurcation with 
8 = 1.00 and a = 1/3 (6 = 30°) 
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Table B.18 


Dependence of y on Re for 


a Round Edged Bifurcation with 


8 = 1.00 and a = 5/9 (0 = 50°) 


is 


Oe 


nq ey - 2. b rho eee 
Le A ‘ sh 


oy 00, fog. 


<oatis in 
ae Gar i Py Dey 


Me aa 


ae 
y : = ae aig a if i 
recipe | : 


1 
an 


105 


Tablet B.19 


Dependence of y on Re for 


a Round Edged Bifurcation with 


8 = 1.00 and a = 7/9 (6 = 70°) 
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Table B.20 
Dependence of y on Re for 
a Round Edged Bifurcation with 
8 = 1.00 and a = 1.00 (6 = 90°) 
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Table B.2]1 
Dependence of y on Re for 


a Round Edged Bifurcation with 


8 = 0780 ‘and a = 1/73 (6 = 30°) 
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Table B.22 
Dependence of y on Re for 


a Round Edged Bifurcation with 


B = 0.80 and o = 5/9 (6 = 50°) 
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Table B.23 


Dependence of y on Re for 
a Round Edged Bifurcation with 
8 = 0.80 and a = 7/9 (6 = 70°) 
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Table B.24 
Dependence of y on Re for 
a Round Edged Bifurcation with 
B = 0.80 and a = 1 (@ = 90°) 
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Table B.25 
Dependence of y on Re for 
a Round Edged Bifurcation with 


8 = 0.60 and a = 1/3 (6 = 30°) 
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Table B.26 
Dependence of y on Re for 
a Round Edged Bifurcation with 
B = 0.60 and a = 5/9 (6 = 50°) 
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Papie’ B27 
Dependence of y on Re for 
a Round Edged Bifurcation with 
8 = 0.60 and a = 7/9 (6 = 70°) 
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Table B.28 
Dependence of y on Re for 


a Round Edged Bifurcation with 


B = 0.60 anda=1 (6 = 90°) 
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Table B.29 
Dependence of y on Re for 
a Round Edged Bifurcation with 
B = 0.40 and a = 1/3 (6 = 30°) 
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Table B.30 
Dependence of y on Re for 


a Round Edged Bifurcation with 


8 = 0.40 and a = 5/9 (0@ = 50°) 
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Table B.31 
Dependence of y on Re for 


a Round Edged Bifurcation with 


8 = 0.40 and a = 7/9 (6 = 70°) 
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Table B.32 
Dependence of y on Re for 


a Round Edged Bifurcation with 


8 = 0.40 anda=1 (6 = 90°) 
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